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Abstract 

Electron-capture and /3-decay of nuclei in the core of massive stars play an 
important role in the stages leading to a type II supernova explosion. Nuclei in 

>, 

CN , the f-p shell are particularly important for these reactions in the post Silicon- 

es 
O 
ON 
O, 

, energy distribution of the Gamow- Teller Giant Resonance (GTGR) for mid- 

r-j ! fp-shell nuclei in terms of a few shape parameters, using data obtained from 



burning stage of a presupernova star. In this paper, we characterise the 



high energy, forward scattering (p,n) and (n,p) reactions. The energy of the 



o 

c/2 ! GTGR centroid Eqt is further generalised as function of nuclear properties 

like mass number, isospin and other shell model properties of the nucleus. 
Since a large fraction of the GT strength lies in the GTGR region, and the 
GTGR is accessible for weak transitions taking place at energies relevant to 
the cores of presupernova and collapsing stars, our results are relevant to the 
study of important e~-capture and /3-decay rates of arbitrary, neutron-rich, f- 
p shell nuclei in stellar cores. Using the observed GTGR and Isobaric Analog 
States (IAS) energy systematics we compare the coupling coefficients in the 
Bohr-Mottelson two particle interaction Hamiltonian for different regions of 
the Isotope Table. 

Ms number PACS number: 97.10.Cv, 23.40.-S, 97.60.Bw 
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I. INTRODUCTION 



It is well known [|l|], 0] that during stellar collapse, the final mass of the homologously 
collapsing core and the strength of the subsequent type II supernova shock is determined 
by the final electron fraction of the core Y e f. The latter, in turn, is influenced by the 
electron captures (and /3-decays) taking place on the nuclei present in the core during it's 
hydrostatic evolution, Kelvin- Helmholtz contraction and dynamic collapse phases. This can 
be seen from the scaling relations JT|, H for the mass of the homologously collapsing core 
and for the shock energy when the core bounces at supernuclear densities:- 

M HC oc Y e 2 f 

and 

E shock ~ ^^{Y ef - Y ei ) ~ Mjg(y e/ - Y ei ) ~ Y e f 3 (Y ef - Y ei ) 

where Mhc and Rue are respectively the mass and radius of the homologously collapsing 
core. The electron-captures and /9-decays also determine the physical conditions in the 
(quasi-) hydrostatically evolving core through their influence on the entropy per nucleon Si. 
This is because a low value of Si implies that (a) the number of free protons is smaller and 
(b) the nuclei are excited generally to low energy states. Since the core loses energy (through 
neutrino emission) mainly by electron captures on free protons, as a consequence of (a), the 
electron fraction Y e remains high leading to a more massive homologously collapsing core 
and subsequently a more energetic shock. Further, since the collapse is essentially adiabatic 
and the nuclei have low excitation energies, the low value of Si ensures that the number 
of drip nucleons is very small, i.e. the nucleons remain inside the nuclei right until the 
core reaches nuclear density. Thus since the number of free nucleons is low, the collapse 
proceeds to higher densities leading to a stronger bounce shock. For all these reasons, to 
make a realistic model of the physical conditions prevailing in the core of a type II supernova 
progenitor, it is necessary to know the weak interaction rates of a number of nuclei at the 
relevant energies. 
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In general, the energy dependence of weak interaction matrix elements (or equivalently, 
the Gamow- Teller (G-T) strength distribution) is unknown for many nuclei of potential im- 
portance in pre-supernova stars and collapsing cores. However for the few nuclei of interest 
for which the experimental data is available, it is possible to map out the energy distribution 
of the GT strength in terms of a few shape parameters. These parameters can be related 
to nuclear properties like the nuclear isospin and the mass number and a corresponding ap- 
proximate distribution can be constructed for any nucleus of astrophysical interest in that 
nuclear shell. We consider here nuclei in the f-p shell since in the post-Silicon-burning phase 
the core of a pre-supernova massive star has a significant abundance of these neutron rich 
nuclei with A > 60. In e~-capture and /3-decay calculations relevant to the pre-supernova 
scenario, it is not particularly important to consider the higher shells (though this is not 
the case during the collapse phase, where at sufficiently high density they may have non- 
negligible abundances). This is because neutron shell blocking [Q] (which starts at around 
74 Ge) substantially decreases the e~-capture rate for these heavier nuclei and their contribu- 
tion to the collective electron capture rates is small compared to that of mid f-p shell nuclei, 
if the latter is present with sufficient abundances. 

The importance of knowing the centroid of the GT distribution (Eqt) li es with the fact 
that it determines the effective energy of the e~-capture and /3-decay reactions from e.g. 
the ground state of the initial nucleus to the excited state of the final nucleus, and this 
along with the electron- Fermi energy determines which nuclei are able to capture electrons 
from, or /3-decay onto the Fermi-sea at a given temperature and density — thus controlling 
the rate at which the abundance of a particular nuclear species would change in the pre- 
supernova core. For example, the following expression holds for the /3~-decay rates of the 
discrete states of the mother nucleus (see e.g. 
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where 

G = ^(2Ji+l)eM-Ei/kT) 

i 

gA and gy are the axial- vector and vector coupling constants respectively and the \M l FGT \ 2 
are the corresponding Fermi and GT matrix elements. E{ is the energy of the i th level of 
the mother having spin J«, and Q- value Qi = Q + Ei. The function f(Qi — E ) is the phase- 
space factor, integrated over electron-energies (in units of m e ) from 1 to {Qi — E ). Eq. (1) 
shows that the distributions of the transition strengths in energy are important in the rate 
calculations. 

The /3-decays can take place either through the Fermi (vector) type interaction or the 
Gamow- Teller (axial-vector) type interaction while the e~-captures on a nucleus in ground 
state involves only the GT interaction. The GT operator G a J2i &itf does not commute 
with the strong spin, isospin dependent forces of the nuclear Hamiltonian, which causes a 
mixing of states in both the spin as well as isospin space. While the effect of mixing in the 
isospin space is small (because of the relatively weaker Coulomb potential term), the mixing 
of states in the spin-space gives rise to a broad distribution of the total GT strength in 
excitation energies known as the Gamow- Teller Giant Resonance (GTGR) which contains a 
large fraction of the total weak-interaction strength sum. In contrast, the Fermi strength is 
concentrated in a narrow region of excitation energies because the Fermi operator commutes 
with all parts of the nuclear Hamiltonian except the Coulomb part. At temperatures and 
densities (T ~ 0.5MeV and p ~ 10 w g/cm 3 ) characteristic of the presupernova core, the 
main contribution to the weak interaction rates is expected to come from the GTGR region. 

In sections |TJ and |TTT] we use the available data on nuclear charge-exchange reactions 
on mid fp-shell nuclei to characterise the GT strength distribution for arbitrary nuclei of 
astrophysical interest. In section |IV| we compare our results with other theoretical methods, 
namely the Ml method used by Klapdor |15j and the method of Fuller, Fowler and Newman 
|T3| . In section we discuss the implications of the observed GT energy systematics vis-a-vis 



the Bohr-Mottelson two body Hamiltonian and in section VI we give our conclusions. 
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II. THE GAMOW-TELLER STRENGTH DISTRIBUTION 



It is well known [[?[], || that the information on charge exchange reactions in nuclei can 
be used to extract the GT strength distributions with respect to the excitation energy 
of the daughter nucleus. The charge-exchange (p,n) and (n,p) reactions are iso-spin and 
spin-dependent over a wide range of projectile energies. The weak interaction strength 
distribution (Bq T or Bq T ) is propotional to the AL = 0, forward scattering cross-section in 
the high projectile energy, low momentum transfer limit. This can be seen by comparing 
the corresponding operators |T2|]:- 



X! V aT {r ip )(Ji.(j p Ti.Tp and ^ V T ( 

i i 

which are similar to: 

i i 

The AL = scattering cross-section for various mid-fp shell target nuclei undergoing (p,n) 
and (n,p) reactions are reported in ref. [16| to f26fl . These cross-sections are extracted 



out of the Odeg spectrum by means of multipole analysis using Distorted Wave Impulse 
Approximation calculations (see ref. Jl2| and references therein). The cross-sections can be 
analysed to obtain the experimental distribution of GT strength in energy to a reasonably 
good accuracy by using known calibration relations (e.g. the relation developed in |14| ) 



and this method has been widely used in the literature (JT6], |]22 



The calculation of the weak-interaction mediated reactions under the astrophysically rel- 
evant conditions requires the knowledge of not only the total strength in a given direction 
( e.g. /3~-decay or e~-capture) but also that of the strength distribution in nuclear exci- 
tation energy. Attempts have been made to obtain these distributions theoretically from 
shell-model calculations. However, the number of shell-model basis states can get very large 
for mid-fp shell nuclei, even at low energies. Hence the direct method of obtaining GTGR 
distribution from shell model calculation using full Ohu basis is computationally very in- 
volved even for a few nuclei. Several attempts to calculate the GT distribution from direct 



shell model calculations using a truncated shell model basis space have also been made, (e.g. 
Aufderheide et al. f|] ) which because of the above reasons, despite substantial computational 
efforts, are as yet approximate. We note that complete fp-shell Otkv Monte Carlo calculations 
have been mentioned recently which show significant quenching (discussed later), and have 
reportedly reproduced the experimentally observed strength in nuclei such as 54 Cr, 54 Fe, 
55 Mn and 56 Fe but not in 58 Ni, the result being sensitive to the interaction Hamiltonian 



assumed (see Koonin and Langanke [TT[] and references therein). In any case the calculation 
of the weak interaction mediated reaction rates for a moderately large number of nuclei 
required in astrophysical situations requires a straightforward and computationally man- 
ageable approach. Here, following an earlier work || we use the framework of a statistical 



approach, as in the Spectral Distribution Theory , [TIJ . 

The sum rule strength for a transition from an initial state \i) to a final state |/) is given 

by 

Sot = E K/i^)| 2 = EMM^/X/Mi) = (*|(0T)W|i) 

/ / 

On the single particle level, if (n^) and (n^) are the fractional occupancies of the neutrons 
in nlj level and protons in nlj' level, then the GT sum rule (e.g. for /3~-decay) is given by 
the expression: 

S%I = 3Z n £ |C^'| 2 (1 - <n£.,»«> (2) 

nljj' 

where C J J = (-1) { J - j)[2(2j + l)(2j' + \)]V 2 W and W is the Racah coefficient. 
The occupation numbers (n^.,), (n^-) etc. in a given shell can be calculated with the 
Spectral Distribution Theory ( see e.g. ref. [[J and ]nj) 

For the distribution of the strength with the nuclear excitation energy we note that the 
according to the Spectral Distribution Theory, for a many particle space of large dimensional- 
ity, the smoothed-out eigenvalue distribution of a (2+1) body Hamiltonian is approximately 
a Gaussian. A skewed Gaussian (called the Edgeworth expansion) of the form: 

B% T = Aq[1 + 7i (x 3 - 3x)/6 + 7 1 2 (x 6 - 15x 4 + 45x 2 - 15)/72] exp(-x 2 /2) 
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where 

x = (E ex - E GT )/(j 

is used here. The parameters 71 (the skewness factor), Eqt (the energy centroid), a (the 
effective half-width) and Aq (normalisation factor) are obtained for each daughter nucleus 
by fitting the above formula to the experimentally obtained strength distribution. 

It is known that if the total strength in the GT + direction is small, the observed value of 
the total Gamow- Teller strength in the GT~ direction is, on the average, approximately 50% 
of the theoretically predicted value of 3(iV — Z) ( for the nuclei studied here, it is observed 
to vary between 45 to 63%). This quenching of GT strength given by the factor Z n in eq. 
(2) above is expected to be due to the excitation of A-isobars at higher excitation energies 
through N^ 1 — A transitions, or the missing GT strength may lie at higher excitations 
(between 30 MeV to 50 MeV) ||. On the other hand in the electron capture direction, 
the total GT strength sum ( i.e. B~Q T (E')dE' ) in the f-p shell has been argued to 
depend on the number of valance protons and the number of neutron holes in the fp shell 
(Koonin and Langanke, 1994). According to these authors the total GT + strength observed 
experimentally in the mid-fp shell nuclei can be fitted by an expression like: 

B(GT+) = 0.0429Z„ ai (20 - N val ). 

As argued by them, in the electron capture direction for nuclei having N > Z, the number 
of possible transitions are limited by the number of available neutron holes and the number 
of protons in the full f-p shell. They also suggest that the neutron-proton correlations 
throughout all of the f-p shell are important in reproducing the observed quenching, and 
mid-fp shell nuclei in this sense behave as though there is only one large shell in which all 
subshell structure has been diluted. The above expression includes the observed quenching 
and can therefore be used to independently fix Z n for the sum-rule for electron capture in 
an equation similar to eq. (2) 

For reasons discussed in section | and the beginning of this section, astrophysical calcula- 
tions require, apart from the total GT + strength itself (such as given above), it's distribution 
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in energy — which is characterised by higher moments of the distribution such as the energy 
centroid E G r, the skewness of the distribution, etc. With the objective of generating the 
Bgt{E) distribution of any arbitrary nucleus (A,Z) of potential astrophysical importance, we 
now make use of the experimentally available data. The Bq T (E) and Bq T (E) distributions 
obtained from analysis of charge exchange experiments are least-square fitted to the Edge- 
worth expression discussed above, under the constraint that the total 'experimental' GT 
strength is reproduced by the area under the fitted expression. The Edgeworth expansion 
gives a good fit in the cases of 51 V, 81 Br, 71 Ga and 60 Ni. In the case of 58 Ni (target nucleus), 
the resonance region between 7.5 and 12 MeV is fitted well by the Edgeworth expression. 
Similarly for 54 Fe, the fit is reasonably good between 6 to 12 MeV. For these nuclei, the 
strength fluctuates rapidly at low excitation energy, suggesting that at these energies single 
particle transitions dominate over collective resonance. Since Spectral Distribution Theory 
is based on a statistical description of a large number of basis states, it is not expected that 
the Edgeworth expansion will reproduce the fluctuating strengths at low energies due to 
single particle transitions. Nevertheless for the G-T collective resonance region (which is 
important because of the bulk of the total strength being accessible at energies relevant to 
the astrophysical situation), the Spectral Distribution Theory approach is useful. 

Finally, since for 56 Fe, 54 Fe, 58 Ni, 60 Ni, only the a(L = 0) (p,n) cross-section is reported 
in ref fll6f, the calibration relation based on the factorised DWIA expression for L = 



differential cross-section as developed in [|1J] was used to convert the reported cross-section 
to BQ T /MeV . Table 1 gives the Edgeworth parameters for each of these fits. A typical fit 
of Bq T per 0.1 MeV interval (for 71 Ga) is shown in Fig. 1. The data used in these fits came 



from refs. || ( 54 < 56 Fe, 58 - 60 Ni), @ ( 81 Br), @ ( 51 V) and @ ( 71 Ga). 



The Edgeworth parameters for Bq T (E) per MeV vs. E x fits are tabulated in Table [TI|. 
Fig. 2 shows the distribution and fit for 59 Co. The distribution is well reproduced for the low 
energy region and the main GTGR peak. For this set of results the data was obtained from 
refs. ( 54 Fe), @ ( 56 Fe, 58 Ni, 55 Mn), @ ( 51 V, 59 Co), @ ( 70 Ge). For 70 Ge, Vetterli et al. 



2~4fl quote the value of Bq T upto 7.6 MeV only, because of the difficulty of extraction of the 
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AL = component from the (p,n) cross- section. The parameters for this nucleus reported 
in Table [Tl] are obtained on the basis of the strength upto 7.6 MeV only. In some cases the 
single-particle transitions may dominate in either direction at low excitation energy, and so 
wherever they are known experimentally, these should be used explicitly in the transition 
rate calculations. 



III. PREDICTED CENTROIDS OF GT STRENGTH DISTRIBUTIONS 

We refer to the centroids of the Bq T (E) distributions obtained by making Edgeworth fits 
to those reported from the (p,n) and (n,p) experiments as the 'experimental centroids'. In 
order to be able to predict such an energy distribution for any relevant f-p shell nucleus, 
we relate the shape-parameters for Bq T distributions to nuclear properties like the nuclear 
ground-state isospin, and the mass number. Here we report the relation developed for the 
energy centroid Eqti which can be applied to arbitrary nuclei of astrophysical interest for 
rate calculations. 



A. GT energy systematics 

The GT ± operator is a space vector and an isovector, and the selection governing this 
transition require that AJ = 0, ±1, An = (no + —>■ + ) and AT = 0, ±1. Since the Fermi 
operator is a pure isovector, the selection rules require that A J = 0, An = and AT = 0, 
i.e. transitions take place only between Isobaric Analog States (IAS). For both Fermi and 
allowed GT transitions AL = 0. It has been argued (see ||, and that under the 



action of resonant GT or Fermi operators, the collective states have centroids located at 
Eq T and Ejas with respect to the appropriate, unperturbed state such that the difference 
Egt- — Eias should depend on the spin-orbit splitting term (~ A^ 1 ^ 3 ), and on the isospin 
dependent Lane potential term (~ (N — Z)/A). The latter gives the energy difference 
between two levels having the same T z , but differ in T by a one unit. The excitation energy 
of the Isobaric Analog State Eias (where the Fermi centroid is located) is obtained in our 
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analysis, either experimentally (where available), or from the following theoretical relation 



developed by Fowler and Woosley and reported in ref |13| for neutron rich nuclei : 



1 72S( Z 1") 

E IAS = AM a - AM C - 0.7824 + ^ -MeV (3) 

R 

with 

R = 1.12A 1/3 + 0.78 

Here AM^ is the mass excess of the initial nucleus, AMq is the mass excess of the final 
nucleus, and R is the radius in fm. This is observed to agree very well with the experimental 
data used in this work (see Fig. 3). Fitting a linear combination of the spin-orbit and isospin 
dependent terms to the experimental centroids, Eq T for the (p,n) reaction data, we get the 
following empirical relation :- 

Egt _ _ ElAS = U.1QA- 1/3 - 76.1(N - Z)/A. (4) 

The correlation between the theoretical E GT - obtained from (4) (together with the experi- 
mental values of Ej^s reported in the references given earlier), and the experimental GT~ 
centroids as referred to above is shown in Fig. 4. It is noted that a similar relation derived 
earlier || (also in 0) using data for nuclei well beyond f-p shell had the corresponding 
coefficients differ substantially from those obtained in eq. (4) for the f-p shell nuclei. The 
reason for this difference is discussed in section |V|. The goodness of the fit is measured in 
terms of the rms deviation from the experimental value and is evaluated to be 0.43 MeV. As 



it is apparent from Table |TTT| and Fig. 4, the GT~ centroids are located roughly between 9 
to 11 MeV. This, together with the expected thermal spread of the electron Fermi distribu- 
tion at temperatures relevant in the presupernova core (typically 5 x 10 9 i^), and a on — 3 
MeV would make a typical error of the order of 0.43 MeV in the predicted centroid quite 
acceptable for astrophysical rate calculations. 

B. GT + energy systematics 

Since in the n-p reaction there is a transition from a parent nuclear ground-state having 
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isospin T = T z to a daughter state having minimum isospin T d = T z + 1 (for g.s. to g.s. 
transition ), there can be no IAS state in the daughter corresponding to the ground state of 
the mother for the electron capture or the GT + direction. Thus, E GT + has to be obtained 
directly (in contrast to E GT - — E IAS in the GT~ direction) and can be expected to depend 
on the spin-orbit splitting term, the Lane potential term and, for odd-A nuclei on the pairing 
energy term as well. 

Electron-capture/e + -emission on a nucleus generally raise the final nucleus to an excited 
state. The single-particle excitation configurations of most of these excited states can be 
constructed from its ground state by breaking either a proton or a neutron pair and raising 
a single particle to an excited level or by raising an unpaired single particle (if available) to 
a excited state. For odd-even/even-odd nuclei, the excited states of the final even-odd/odd- 
even nuclei (which are connected to ground state of the parent nucleus by GT + transitions), 
can be generated from the g.s. of the final nucleus by breaking a particle pair in lower 
energy level and using it to create an unpaired particle in a higher energy level or to pair 
off a previously unpaired particle in a higher energy level and simultaneously raise another 
single particle to the excited state. For an odd-odd / even-even initial nucleus, however, the 
excited states in the final even-even/ odd-odd nucleus will be generated by breaking and 
making a particle pair. This implies that the energy required to break a pair in the final 
odd A nucleus must also be included in any expression for GT + energy systematics. 

Fitting the experimental E GT centroids (obtained from the Edgeworth weighted exper- 
imental B GT (E) distributions) to a linear combination of these terms, we arrive at the 
following empirical relation for E GT + centroid with respect to the ground state of the final 
nucleus: 

E GT+ = 13.1(L4- 1/3 - 11.28(iV - Z)/A + 12A~ 1/2 5 Aodd . (5) 

The plot of E GT + from (5) vs. the experimental centroid is shown in Fig. 5. The goodness 
of fit in Fig. 5 is also measured in terms of the r.m.s. deviation from the experimental 
value, and has the value 0.31 MeV. The nuclei used in both figures are stated in the order of 
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increasing Eqt (experimental). We note that Koonin and Langanke (1994) suggest that the 
GT-resonance appears in (n,p) spectra at systematically higher energies for odd-Z targets 
than for even-Z targets, although this is based on data involving target nuclei with even 
number of neutrons. On the basis of single particle excitation diagrams for the final nucleus, 
conforming to G-T selection rules for a transition from the initial to the final nucleus, we 
find that the pairing energy dependence in terms of odd-A in eq. (5) (the last term) would 
be present even in odd-N, even-Z nuclei. 

The average difference of experimental and predicted values of 0.4 MeV in Eq T and 
0.3 MeV in the case of Eq T can be mitigated by the spread in the electron Fermi-Dirac 
distribution due to the 0.5 MeV temperature and also due to the fact that in a given 
region of the star at a particular stage of evolution there are several nuclear species present 
in the stellar material simultaneously. The total transition rate of this large admixture 
of nuclei is likely to be dominated in a given stage and region by a few nuclear species, 
with relatively large transition rates and abundances. These would necessarily have large 
energies of transition and would in any case consume a substantial fraction of the strength 
sum rule. Thus for astrophysical purposes, the level of accuracy of the predictions for the 
GT 1 * 1 centroids obtained here are expected to be adequate. 

Data for 45 Sc and 48 Ti are also available in references |2lJ and p6 [. However, we have 
not used these two nuclei in obtaining eq.(5) and Fig. 5, partly because these nuclei are 
close to the beginning of the fp-shell where the general assumptions of the statistical nature 
of the Spectral Distribution Theory may not be valid. For this reason, the calculation of 
the Edgeworth weighted centroid of the experimental Bqt+ {E) distribution gives inaccurate 
results because of the low energy strength. Indeed, if these nuclei were to be included in 
Fig. 3, the goodness of fit would have decreased to an r.m.s. value of 1.37 MeV. Similarly, 
we had to ignore the data for 90 Zr [^5[ because the g 9 / 2 shell is occupied in this nucleus. 
Effectively, it can be said that eq. (5) is valid for mid-fp-shell nuclei with N > Z. 
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IV. COMPARISON WITH OTHER THEORETICAL METHODS 

While a detailed shell model calculation for these nuclei requires a very large number of 
basis states and would be difficult to carry out for each of these nuclei, it is conceivable that 
a reasonably good estimate of the GT excitation energy centroids by using either the Ml 
excitation method (for the GT + excitation) or the FFN method (for the GT~ excitation). 
These methods are described in detail in ref. ||13[| . The predictions of these methods are 
compared with the results obtained in this paper in the following paragraphs and in Table [Hi 
and Table [TV]. 

In the FFN method, the GT transition in GT + (GT~) direction, which can be either 
the spin-flip (sf) or the no-spin-flip (nsf) transition, is reproduced by taking into account all 
possible excitation of a p(n) in the parent nucleus to a sf or nsf orbital (for which AL = 0) 
and then converting the p(n) to n(p) via the t + (t~) operator. The GT excitation energy for 
that transition is given by the sum of the excitation energies with respect to the ground state 
of the daughter nucleus together with the Lane potential energy and if required, the pairing 
energy term. The GT centroid energy with respect to the ground state of the daughter 
nucleus is given by the weighted average of all these sf and nsf transitions. For example, in 
Fig. 6, the ground state of (g.s.) 51 V is connected to the g.s. of 51 Cr by a nsf-r~ transition, 
and to an excited state by a sf-T~ transition. For the former transition, the GT excitation 
energy (measured, as always, with respect to the ground state of 51 Cr) equals the Lane 
potential ( 50(iV — Z)/A), which is evaluated as 4.90 MeV. The square of the GT matrix 
element for this transition is 6.42. For the latter transition, the GT excitation energy is the 
sum of single-particle excitation energies (6.66 MeV), the Lane potential (4.90 MeV) and 
pairing energy (1.68 MeV), while the square of the transition matrix element is 13.71. Here, 
the sf transition has the maximum contribution to the total strength and this is observed 
to be the norm in the case of other nuclei considered in ref. (T^j as well. Thus, from this 
single particle excitation method, the GT~ centroid can be expected to lie at 10.58 MeV. 
The experimental centroid lies at 11.75 MeV. However, as this method does not take into 
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account the collective particle-hole excitations || of the nucleus (where a good fraction of 
the total strength can be expected to lie), it should not be expected to yield the experimental 
Eqt precisely. This is a purely theoretical method, and the excitation energies are calculated 



from single-particle shell model energy levels [[!(J. For the single particle levels of neutron 
we have used the sequenece and values labelled as "Seeger" in ref [[!tj. Although for the 
protons, the "Seeger" sequence lists 2p 3 / 2 higher than H5/2, since the observed ground state 
spins of 63 Cu, 65 Cu etc (where there is one extra proton beyond the H7/2 level) contradict 
this sequencing (i.e. the g.s. J n for these nuclei is (3/2)~ and not (5/2)~), following an 
earlier work || we adopt the same values for the proton single particle energy differences as 
in the case of neutrons in this work unless specifically mentioned otherwise. 

Now, it is seen from the selection rules that most of the GT Strength for nuclei with 
A ~ 60 lie in the spin-flip transitions. Experimentally too, it is observed that the sf transition 
strength is more concentrated than the nsf transition strength. For the sf transitions, the 
action of the GT operator is equivalent to a s-f transition followed by an isospin-flip (Ml 
transition from the mother to the daughter). If we consider a T < to T > transition, to find 
E GT + in the daughter by the Ml method, we need to consider the Ml-r~ transition from 
the daughter to the mother nucleus. The Ml method for finding the GT + -centroid in the 
neutron rich daughter nucleus (see e.g. |l3j ) is based on the observation that most of the GT 
sf- configuration strength is concentrated near the Anti-Isobaric Analog State (AIAS) of the 
Ml giant resonance — the so called Ml-AIAS state. Now, the predominant configuration in 
the Ml-AIAS state in the daughter is the spin-flip configuration generated from a T < mother 
ground state by transforming a proton in the j = 1 + 1/2 state into an empty neutron in 
the j = I — 1/2 level. The spin-flip configuration in the daughter is therefore the starting 
point and it is this configuration whose excitation with respect to the g.s. would yield the 
sought for GT-centroid. The Ml method generates, in the T < mother, the T > analog of 
this sf- configuration by the application of the T + = Y^t=\ T t (isospin-raising) operator . In 
general, there may be more than one AIAS of the sf-configuration which are orthogonal to 
each other. These AIAS states contain most of the Ml excitation configuration strength 
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and their energy can be calculated from the single-particle shell model, taking into account 
the particle-hole repulsion energy and the pairing energy wherever relevant. The T > Analog 
state is seperated upwards from the AIAS states by the Lane potential. Once the energy 
of this T > Analog state in the mother ( of the sf- configuration in the daughter ) is known, 
subtraction of the energy of the Analog state corresponding to the g.s. of the daughter ('the 
first analog') would yield the approximate excitation energy of the sf-configuration in the T > 
daughter with respect to its g.s., which is the required GT + centroid in the daughter. This 
assumes an argument similar to the Brink hypothesis, as discussed later. The first Analog 
state in the mother is often known experimentally, and where available, this is an advantage 
in the Ml-method. 

As an illustration of this method, we consider the case of GT + centroid in 54 Mn (for 
the 54 Fe(e~, i/ e ) 54 Mn reaction). Fig. 7a shows the g.s. configuration of the mother nucleus 
54 Fe. The spin-flip excitation of 54 Mn are shown in Figs. 7b and alongside it are shown the 
excited states of 54 Fe which are obtained by operating r~ on the spin-flip state of 54 Mn. 
Thus the GT excitation corresponding to the sf excitation in Fig. 7b (in 54 Mn) will be a 
superposition of two shown basis configurations. Now, the AIAS corresponding to these two 
configurations can be constructed and its position in energy is determined as the weighted 
sum of the single particle excitation energies and the Lane potential (3.70 MeV for this case). 
The weighting factors obtained from Ml method for the 54 Fe excited state configurations are 



y3/4 and y 1/4 respectively (see Fig. 7) and the excitation energies for each configuration 
in Fig. 7b are noted below them. The position of the GT excitation corresponding to the 
s-f configuration of 54 Mn will then be given by the difference in the energy of the obtained 
IAS state and the IAS state in 54 Fe corresponding to the g.s. of 54 Mn (the so called first 
analog state — located at an excitation energy of 8.7 MeV — see below). 

A few examples of the Ml method in T < to T > transitions are given in Table [IV]. For all 
the nuclei quoted in Table [IV], the required Ejas in the mother is known from experiment. 
For the mother nuclei quoted in Table [ITI| the Ejas of the daughter ground states were 



obtained from f29| ( 56 Fe) and |28| ( 59 Co) and are 11.51 and 9.55 MeV respectively. For 54 Fe 
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we note that ref. [[27| gives only the Isobaric Analog States corresponding to a few excited 
levels of 54 Mn. The difference of the IAS energies and the corresponding levels in 54 Mn are 
in the range of 8.59 to 8.79 MeV. Using the assumption (similar to the Brink hypothesis in 
electromagnetic transitions — see that the energies in the resonance states scale in the 
same way as those of the discrete states in the daughter, we adopt the mean value of 8.7 
MeV as the Eias m 54 Fe corresponding to the g.s. of 54 Mn. The numbers quoted in the 
last column of Table |IV| are for IAS states giving GT energy closest to the experimentally 
observed GT-centroid . 

For the nuclei under consideration, the experimental inputs coupled with the given the- 
oretical prescriptions (Ml and FFN method) give numbers which can be compared with the 
empirical relations obtained here. Tables |T| and [TV] show that the empirically fitted values 



of Eqt+ and Eqt- as reported in section 111 are somewhat closer to the experimental cen- 



troids than the numbers obtained from the Ml and FFN methods. The results of sections |T( 
and [Til] can therefore be used to calculate the weak transition properties of f-p shell nuclei 
of interest to astrophysics. 

V. GT COLLECTIVE STATES FROM DISPERSION RELATIONS WITH 
BOHR-MOTTELSON HAMILTONIAN 

Attempts have been made in the past to obtain a simple mass formula for the GT~ 
energy systematics using the single particle level structure and two particle Bohr Mottelson 
type interaction Hamiltonians. The GTGR and the IAS states can be considered to be 
collective excitations excited by the spin-dependent and charge-exchange Y*t=\ T-^ 1 an d 
J2t=i T- operators. More specifically, in the GT~ direction, the energy difference Eqt — Ejas 
can be written as ( see e.g. Suzuki pi| , 1982): 



M 




7r) (tt\ 


Eiil t+H Y,f=i tI 




(7T 




7r) (7T 







Eqt- — Ejas 

The following form of the two particle interaction Hamiltonian (Bohr and Mottelson, 1981) 
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for example has been used here by various authors (e.g. ref. ||, |52] 



A 1 1 1 

H = - ^UiOi + -Kr Ti.Tj + -K a Oi.Oj + -K aT Y( T i- r j){ a i- a j) 
i=l i^j i^j - i^j 

where « T and the coupling constants (divided by the mass number) for the 

spin-spin, isospin-isospin and the spin-isospin interaction terms and & is the single particle 
spin-orbit coupling constant. Using the simplification brought about by considering the 
case of very neutron rich nuclei where the Tamm-Dancoff approximation would be valid 
and further considering only nulclei for which the protons occupy the major shells and the 
neutrons fill the j — 1 + 1/2 subshell, e.g. in the cases of 90 Zr and 48 Ca, Suzuki |H| has used 
the experimentally known values of Eqt- ~ Eias to determine the difference in the coupling 
constants ( Table 6). 

Another approach (due to Gaarde ||) to the GTGR energy systematics is from the field 
description of the coherent particle-hole excitations which lead to the GT resonance. In the 
field description, the coherent state is generated by an oscillation of average field in the spin- 
isospin state which is propotional to err. The self-consistency condition on this oscillating 
potential leads to a dispersion relation for the energy of the GT collective state. Further, if 
it is assumed that the strength is clustered mainly around two regions, one corresponding 
to the no-spin-flip transitions with energy e« and the other to the spin-flip transitions with 
energy 6j + A; s then, the energy dispersion relation equivalent to an RPA equation for a 
seperable force can be written as: 

2(N-Z)(l-f) + 2{N-Z)f = -1 
e?: — e Q + Aj s — e n UT 

where / is the fraction of the total strength lying in the spin-flip region, and Aj s is the 

spin-orbit splitting energy. The quantity / is evaluated for each of the nuclei (for which 

(p,n) data exists) from the single particle transition configurations constructed using the 

FFN method and is tabulated in Table 0. A similar expression may be developed for the 

IAS transition which has the form: 

eiAs ~ djAs = 2k t (N - Z) (7) 
17 



Using the experimentally known values of Eias an d eq. (7), we obtain the value of k t 
as 49.04/A Keeping this value of k t fixed, the 'experimental' values of Eq T and Ejas 
obtained from Table [Til] and Fig. 3 are fitted by a two parameter least square fit to the 
quantity (ear — U,gt) — (cias ~ e i,iAs) , where the first term is evaluated from the quadratic 
roots of eq. (6). This gives the value of K aT as 11.16/A As is shown in Table [VT|, these 
values differ from those obtained for 208 Pb , or even for 90 Zr region, indicating that in order 
to explain both the IAS and GT energetics, the simple Bohr-Mottelson Hamiltonian requires 
coupling coefficients which are substantially different for different shells. At the same time 
it is not surprising that the relation (4) developed for f-p shell nuclei which can be shown 



to depend on the coupling constants (as in the approach of Suzuki f31fl , |32|), would differ 



substantially from similar relations || developed for a different region of the isotope table. 

VI. CONCLUSION 

Since the Spectral Distribution Theory is essentially a statistical theory requiring the ex- 
istence of a large number of basis states, this approach is expected to be valid only for 
nuclei which are far removed from the closed shell and sub-shell configurations. Using the 
framework of the Spectral Distribution Theory, we have used the experimental data on 
charge exchange reactions on fp-shell nuclei to obtain the energy centroid of the collective 
Gamow- Teller resonances for arbitrary nuclei in this shell. These quantities, together with 



the prediction of total GT strength ( such as given in [11| ) are useful in the prediction of 
reaction rates mediated by weak-interactions in the cores of massive pre-supernova stars. 
We have in eqs. (4) and (5) a dependence of the GT energy centroids in the f-p shell on the 
spin orbit interaction, the Lane potential and the pairing energy. Further, with a view to 
comparing the implications of the experimental results for fp-shell nuclei with earlier work 



on heavier nuclei ( see ||, || and [|31[] ) we have used the experimental values of the GT - 
centroids and the IAS energies to extract the relevant coupling coefficients in the Bohr- 
Mottelson Hamiltonian, for the f-p shell. The fact that they differ from those of the earlier 
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studies for heavier nuclei such as in the region of 208 Pb indicates that apart from the 1/A 
dependence, there may be an intrinsic shell dependence of these coupling constants when a 
simple model like the one due to Bohr-Mottelson has to be used. 

The results of these calculations are being used to predict the neutrino energy spectrum 
during the collapse phase upto the point of z/ e -trapping. In the event of a sufficiently nearby 
supernova explosion, these can be compared with the neutrino spectroscopy results obtain- 
able by future neutrino detectors (e.g. SNO, ICARUS) — thereby revealing clues to the 
important thermodynamic and nuclear conditions of the presupernova core. A preliminary 



discussion of this was reported in ref. [34 



We thank K. Kar and S. Sarkar for the initial impetus to work on this problem. This 
research formed a part of the 8th Five Year Plan (8P — 45) at Tata Institute. 
Note added in the proofs: 

After this work was accepted for publication, a paper on experimental G-T strength 
distributions in 60 Ni, 62 Ni and 64 Ni by Williams et al (Phys. Rev. C 51, 1144 (1995)) came 
to our notice. Using their experimental data, we find that the strength weighted energy 
centroids for these nuclei are at 2.59, 2.56 and 2.8 MeV respectively. This is to be compared 
with our predictions (from eq. (5)) of 2.59, 2.21 and 1.865 MeV. While the predictions of 
eq. (5) (based on the (n,p) data of the earlier seven nuclei) are reasonably good for 60 Ni and 
62 Ni, it does not work well in the case of 64 Ni, possibly because in the latter nucleus much 
of the strength is concentrated in the ground state to ground state transition. 



19 



REFERENCES 

[1] G. E. Brown, H. A. Bethe and G. Baym, Nucl. Phys. A375, 481 (1982) 

[2] H. A. Bethe, G. E. Brown, J. Applegate, J. M. Lattimer, Nucl. Phys. A324, 487 (1979) 

[3] G. Fuller, Astrophys. J. 252, 741 (1982) 

[4] M. B. Aufderheide, S. Bloom, D. A. Resler, G. J. Mathews, Phys. Rev. C 47, 2961 
(1993) 

[5] K. Kar, A. Ray and S. Sarkar, Astrophys. J. 434, 662 (1994) 

[6] K. Nakayama , A. Pio Galled and F. Krmpotic, Phys. Lett. 114B, 217 (1982) 

[7] D. J. Horen et al, Phys. Lett. 99B, 383 (1981) 

[8] C. Gaarde, Nucl. Phys. A396, 127 (1983) 

[9] J. B. French and V. K. B. Kota, Ann. Rev. Nucl. Particle Sci.32, 35 (1982) 
[10] V. K. B. Kota and K. Kar. Pramana 32, 647 (1989) 
[11] S. E. Koonin and K. Langanke, Phys. Lett. 326B, 5 (1994) 
[12] T. N. Taddeucci et al, Nucl. Phys. A469, 125 (1987) 

[13] G. M. Fuller, W. A. Fowler, M. J. Newman, Astrophys. J. 252, 715 (1982) 
[14] T. N. Taddeucci et al, Phys. Rev. C 25, 1094 (1981) 

[15] H. V. Klapdor, (1976), CERN-Report 76-13, Proceedings Third International Confer- 
ence on Nuclei Far from Stability. 

[16] J. Rapaport et al, Nucl. Phys. A410, 371 (1983) 

[17] M. C. Vetterli et al, Phys. Rev. C 40, 559 (1989) 

[18] D. Krofcheck et al, Phys. Lett. 189B, 299 (1987) 

20 



[19] J. Rapaport et al, Nucl. Phys. A427, 332 (1984) 
[20] D. Krofcheck et al, Phys. Rev. Lett. 55, 1051 (1985) 

[21] W. P. Alford, A. Celler, B. A. Brown, R. Abegg, K. Ferguson, R. Helmer, K. P. Jackson, 
S. Long, K. Raywood and S. Yen, Nucl. Phys. A531, 97 (1991) 

[22] S. El-Kateb et al, Phys. Rev. C 49, 3128 (1994) 

[23] W. P. Alford et al, Phys. Rev. C 48, 2818 (1993) 

[24] M. C. Vetterli et al, Phys. Rev. C 45, 997 (1992) 

[25] K. J. Raywood et al, Phys. Rev. C 41, 2836 (1990) 

[26] W. P. Alford et al, Nucl. Phys. A514, 49 (1990) 

[27] W. Gongqing, Z. Jaibi and Z. Jingen, Nucl. Data Sheets, 50, 255 (1987) 

[28] P. Andersson, L. P. Ekstrom and J. Lyttkens, Nucl. Data Sheets, 69, 733 (1993) 

[29] H. Junde, H. Dailing, Z. Chunmei, H. Xiaoling, H. Baohua and W. Yaodong, Nucl. Data 
Sheets, 51, 1 (1987) 

[30] M. Hillman and J. R. Grover, Phys. Rev. 185, 1307 (1969) 

[31] T. Suzuki, Nucl. Phys. A379, 110 (1982) 

[32] T. Suzuki, Phys. Letts. B104, 92 (1981) 

[33] A. Bohr and B. R. Mottelson, Phys. Lett. 100B, 10 (1981) 

[34] F. K. Sutaria and A. Ray, Proceedings of 'Nuclei in Cosmos', 3rd International sym- 
posium on Nuclear Astrophysics, Editors: M. Busso, R. Gallino and C. Raiteri, AIP, 
(1994), in press. 



21 



FIGURES 

FIG. 1. The best fit Edgeworth distribution superposed on the B GT distribution obtained from 
p-n reaction on 71 Ga at 120 MeV. The excitation energy E x is with respect to the ground state of 
71 Ge. 

FIG. 2. The best fit Edgeworth distribution superposed on the B GT distribution obtained 
from n-p reaction on 59 Co at 198 MeV. 

FIG. 3. Correlation between Ejas as calculated from eq. (3) and the experimental positions 
of the IAS for fp-shell nuclei from (p,n) reaction data. Nuclei indicated in the figure are ordered 
with increasing values of experimental Ejas- 

FIG. 4. E GT - as obtained from eq. (4) plotted against the 'experimental' E GT - centroids. 
Nuclei indicated are given in order of increasing value of experimental E GT - . 

FIG. 5. E GT + as obtained from eq. (5) plotted against the 'experimental' E GT + centroids. 
Nuclei are given in order of increasing value of experimental E GT + . 

FIG. 6. Schematic diagram of single particle excitations illustrating the application of FFN 
method to determination of GT~ strength centroid with respect to the g.s. of 51 Cr in /?~-decay of 
51 V. The level spacings are not drawn to scale. 

FIG. 7. The Ml method for 54 Fe(e~,i/ e ) 54 Mn. Fig.7a shows the g.s. configuration of 54 Fe. In 
Fig. 7b, the diagram on the left represents a spin-flip excitation in in 54 Mn generated from 54 Fe 
g.s. by transforming a l/ 7/2 proton into a l/ 5/2 neutron. The IAS for the 54 Fe configuration is 
a j^/3/4, ^1/4} superposition of the shown basis states on the right; the AIAS state is similarly 
^/3/4}. The AIAS is seperated from the IAS only by the Lane potential. 



22 



TABLES 

TABLE I. Edgeworth parameters for Bq T energy distribution obtained from p-n reactions. 

Nucleus A (B GT /MeV) E gt ( MeV ) a ( MeV ) 7i 

51 V 1.2 10.7 4.2 -1.00 

81 Br 2.0 10.2 3.6 -1.19 

71 Ga 2.0 9.6 3.0 -1.25 

58 Ni 1.0 9.4 1.8 0.11 

60 Ni 0.8 9.0 2.6 1.05 

54 Fe 1.2 8.9 1.8 -0.50 

56 Fe 1.0 9.0 3.5 -0.67 
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TABLE II. Edgeworth parameters for Bq T energy distribution obtained from n-p reactions data. 



Nucleus A (5J T /MeV) E^ T MeV a (MeV) 71 

54 Fe 0.6 2.8 2.1 1.01 

51 V 0.3 3.8 1.3 -1.21 

59 Co 0.6 4.3 1.2 0.43 

70 Ge 0.1 2.1 2.1 -1.53 

56 Fe 0.6 2.3 1.8 0.97 

55 Mn 0.4 4.4 1.9 0.53 

58 Ni 1.0 3.3 1.5 0.59 
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TABLE III. Experimental values of the energy-centroids in MeV for GT excitation in the 
mother nucleus and values obtained from the fit eq.(4) and from the FFN Method. 



Nucleus 


E GT - (MeV) 


E GT - (MeV) 


E GT - (MeV) 




(Expt.) 


(Eq. 4) 


(FFN) 


51y 


10.7 


11.1 


10.6 


81 Br 


10.2 


9.6 


12.0 


71 Ga 


9.6 


9.9 


13.6 


54p e 


8.9 


8.9 


8.7 


58 Ni 


9.4 


9.0 


6.5 


56 Fe 


9.0 


9.7 


12.6 


60 Ni 


9.0 


8.7 


4.9 
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TABLE IV. Experimental values of the energy-centroids for GT + excitation in the 
mother nucleus and the values obtained from eq. (5) and from the Ml method. 



Nucleus 


E GT + (MeV) 
(Expt.) 


E GT + (MeV) 
(Eq. 5) 


Ml config. a (daughter) 


E GT + (MeV) 
(Ml) fe 


54 Fe 

56 Fe 

59 Co 


2.8 
2.3 
4.3 


3.05 
2.62 
3.97 


^/7/2)'( Zy /5/2' P3/ 2 ) 
( 7r / 7/2 )'^/5/2' l/ p3 /2 ) 


5.17(4.07) 
4.44 (4.60) 
6.29(5.34) 



(a) The single particle shell model configurations []3(J for the outermost levels of 
the daughter in the Ml excited state (Seeger neutron levels for both protons and 
neutrons). The levels are given in order of increasing energy, starting with the lowest 
energy, partially-filled obital. 

(b) The numbers quoted within parenthesis are the E GT values from Ml method when 
the seperate 'Seeger' ]30| levels are used for neutrons and protons. 
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TABLE V. The distribution of B GT strength in the spin-flip (sf) and no-spin-flip (nsf) tran- 
sitions. Column 4 gives the fractional strength of the sf transitions. 

Nucleus Bq T (nsf) Bq T (sf) f 

51 V 6.42 13.71 0.68 

81 Br 10.13 14.47 0.59 

71 Ga 6.60 15.04 0.69 

54 Fe 2.57 13.71 0.84 

58 Ni 3.33 16.37 0.83 

56 Fe 5.90 16.37 0.73 

60 Ni 0.0 19.04 1.00 



27 



TABLE VI. Coupling constants k t , k ot for Bohr-Mottelson two particle interaction. 



Authors 



Nuclei 



A.K n 



Gaarde @ 
Suzuki |32| 
Suzuki 
Nakayama 
et al. § 
This work 



208p b 
90 Zr 
208p b 

'Global' fit 

92 Zr tQ 208p b 

f-p shell 



28 
32.5 a 
32.5 a 
32.5 

49.0 



23 

25.9 - 28.6 
28.1 -29.0 a 
23.25 

11.2 



5 

3.92 
3.95 
9.25 

37.8 



(a) Bohr-Mottelson estimate, quoted in |32| . 
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